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Abstract 


In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning 
tunneling microscope (CryoSTM) is proposed. TIA with the tip-sample component 
in CryoSTM is called as CryoSTM-TIA. With transimpedance gain of 1 GQ, the 
bandwidth of the CryoSTM-TIA is larger than 200 kHz. The distinctive feature of 
the proposed CryoSTM-TIA is that its pre-amplifier is made of a single cryogenic high 
electron mobility transistor (HEMT), so the apparatus equivalent input noise current 
power spectral density at 100 kHz is lower than 6 (fA)?/Hz. In addition, “bias-cooling 
method” can be used to in-situ control the density of the frozen DX- centers in the 
HEMT doping area, changing its structure to reduce the device noises. With this 
apparatus, fast scanning tunneling spectra measurements with high-energy-resolution 
are capable to be performed. And, it is capable to measure scanning tunneling shot 
noise spectra (STSNS) at the atomic scale for various quantum systems, even if the 
shot noise is very low. It provides a powerful tool to investigate novel quantum states 
by measuring STSNS, such as detecting the existence of Majorana bound states in the 
topological quantum systems. 


1 Introduction 


High performance transimpedance amplifier (TIA) for cryogenic scanning tunneling mi- 
croscope (CryoSTM) is a key element BP]. Recently, in Ref.{3], it was proposed a 
design of TIA for CryoSTM. TIA with the tip-sample component in CryoSTM is called 
as CryoSTM-TIA. In that design, the CryoSTM-TIA has a transimpedance gain of 1 GQ, 
a bandwidth of more than 300 kHz, and its equivalent input noise current power spectral 
density (PSD) that characterizes the circuit inherent noise is 21 (fA)?/Hz at 100 kHz. The 
CryoSTM-TIA is capable of measuring the scanning tunneling shot noise spectra (STSNS) 
of various quantum systems at atomic resolution, even if the measured noise current PSD 
is as small as a few (fA)?/Hz. In Ref.[3], it was illustrated how the CryoSTM-TIA is 
used to measure STSNS of quantum systems to investigate novel quantum states, such 


cryoliang@qq.com 


as detecting Majorana bound states (MBSs) in iron-based superconductors [4]. In the 
CryoSTM-TIA, cascade the pre-amplifier (Pre-Amp) and the post-amplifier (Post-Amp) 
to form an operational amplifier (OPA) denoted as Macro-OPA. Pre-Amp in Macro-Amp 
is a differential amplifier made of a pair of high electron mobility transistors (HEMTs), 
where the HEMT is a kind of cryogenic GaAs MESFET [5]{6]. However, it is quite difficult 
to select cryogenic GaAs MESFETs with identical performances for pairing [2]. 

In this work, a design for the CryoSTM-TIA with a single HEMT in Pre-Amp is 
proposed. In the CryoSTM-TIA, cascade Pre-Amp and Post-Amp to form an inverting- 
amplifier (Inv-Amp), which is not OPA. The CryoSTM-TIA still has transimpedance gain 
of 1GQ and bandwidth more than 200 kHz. However, its equivalent input noise current 
PSD is only 5.3 (fA)?/Hz at 100 kHz, i.e. 1/4 of that in Ref.[B], since a noisy HEMT is 
reduced. Therefore, the lower tunneling shot noise of quantum systems can be measured 
at the atomic scale with this apparatus. Furthermore, only a single HEMT is used in 
Pre-Amp, avoiding the difficulty of matching two identical HEMTs. And, “bias-cooling 
method” can be used to in-situ reduce the inherent noise of the HEMT. With these 
advantages, this apparatus will be a powerful tool to investigate novel quantum states in 
various quantum systems, in the wider range of applications than that in Ref.|3], such as 
high-T, superconductors [7 [8], topological superconductors [4], MBSs [4 {9} [12], 
Andreev reflection [13} [4], and Kondo effect [15], etc. 


2 Circuit of the proposed CryoSTM-TIA 


The circuit of the proposed CryoSTM-TIA is shown in Figfi| It consists of several com- 
ponents: the single HEMT amplifier part of Pre-Amp shown in dashed box (al) of Figf1| 
the power supply of Pre-Amp in dashed box (a2), Post-Amp shown in dashed box (b), 
the compensated feedback network shown in dashed box (c), and the signal source circuit 
shown in dashed box (d). The two stage amplifier made of Pre-Amp and Post-Amp is 
called as inverting-amplifier (Inv-Amp). Inv-Amp is connected with the feedback network 
to form TIA. The input G of TIA is connected to the signal source circuit, to form the 
CryoSTM-TIA. The components placed in the cryogenic zone are shown in the dotted 
box. The parameters of all components of the CryoSTM-TIA circuit are listed in Table[I] 


2.1 Design of Pre-Amp 


As shown in Fig{i] only a single cryogenic HEMT is used in Pre-Amp. The HEMT in 
this work is CNRS-HEMT (denoted as HEMT) [5] [6] developed by CNRS/LPN in France 
with the excellent cryogenic and noise characteristics. It is capable of operating under 0.5 
K with only 0.1 mW on the ideal operating point of “Vas = 100 mV and Ias = 1 mA”. Its 
parameters are listed in Table j! where e?, is its equivalent input noise voltage PSD and 
is is its equivalent input noise current PSD [6]. 

The single HEMT amplifier part of Pre-Amp is shown in dashed box (al) in Fig/1] 
The source of HEMT is grounded by a resistor Rs and a small variable resistor Rg. 
Rs = R; + Rsı, and a capacitor Cs of 11 uF is in parallel with Rg. Drain O, of HEMT 
is connected to the load resistor R,. The resistor Rp is grounded through the variable 
resistor Rpi. Ry = Ry + Ry. The other end of Ry, denoted as O2, is connected to the 
resistor R2, and L as the other end of Rə is connected to P (the other end of Rı) through 
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Figure 1: Circuit of the proposed CryoSTM-TIA. Single HEMT amplifier part of Pre-Amp 
is shown in dashed box (al), power supply of Pre-Amp in dashed box (a2), Post-Amp in 
dashed box (b), compensated feedback network in dashed box (c), and the signal source 
circuit in dashed box (d). The components placed in the cryogenic zone are shown in 


the dotted box. The parameters of all components of CryoSTM-TIA circuit are listed in 
Table 


a small variable resistor Rp. Ri = Ro = Ry. P is connected to the power supply of 
Pre-Amp. P and L are grounded by the large capacitors Cy and Co respectively. In the 
power supply, Rp = Ry + Ry, and Ry is a variable resistor. C'a is the input capacitance 
of Pre-Amp. HEMT, R,, Rn, Ry, and R, are placed in the cryogenic zone, shown in the 
dotted box. C{ is the capacitance of the cable that connects Gate G of HEMT and the tip 
in the CryoSTM. HEMT is as close to the tip as possible to reduce Cy, which can be less 
than 0.5 pF [I6]. The operating point adjustment for Pre-Amp is shown in Sect {4 The 
Pre-Amp parameters are shown in Table [i] 

For Pre-Amp, when AC voltage V is applied to input G, the AC voltage difference 
between O; and Oz is denoted as Vics Ayp = Vop / V is the voltage gain of Pre-Amp. Since 
the gain-bandwidth-product of the CNRS-HEMT is gm/[27(Cgs + Cga)] ~ 1 GHz, the 
bandwidth of Pre-Amp is more than 30 MHz. In max{gm/(27Cx), 1/(27RgCg)} << f <3 
MHz (i.e. 30 Hz << f < 3 MHz with the parameters listed in Table [1}, 


AvP © —9m Ra, (2.1) 


where Rg = RL/(1 + gaRı). ie., Ayp ~ —20 in 30 Hz<< f < 3 MHz. The input 
capacitance of Pre-Amp is 
Ca = Cy + (1 — Avp)Cga- (2.2) 


Table 1: Parameters of all components of CryoSTM-TIA circuit 


HEMT (CNRS-HEMT) 


Gate-source resistance Ra >10 TQ 
Transconductance gm 40 mS 
Channel conductance gq 1 ms 
Gate-source capacitance C'gs 5 pF 
Gate—drain capacitance Cga 1 pF 
Drain-source voltage Vas 100 mV 
Drain-source current Ias 1mA 
Z 
Equivalent input noise voltge PSD ex eee ee ai 
2 
Equivalent input noise current PSD iZ ara NG? i 
Pre-Amp 
Rs + Rsı 50+5Q Cs 11 uF 
Ry + Rri 15045 Q Rı, Re 1 kQ 
Rp 5+5 Q 
Ci 22 uF C2 11 uF 
PVR LM4050-10 Ry 2 kQ 
Ri + Ra 4400 + 10 2 Vep +15 V 
THS4021 as Rear-OPA 
Qao 94 dB fo 16 kHz 
Ca 1.5 pF Ra 1 MQ 
VRp +5 V Van -5 V 


2.25 (nV)? /Hz in f > 10 KHz 
4 (pA)? /Hz in f > 10 KHz 


Equivalent input noise voltge PSD e2 
Equivalent input noise current PSD i2 


Post-Amp 
Rr 200 kQ Cu, C; 150 pF 
Re1,Rea2 100 Q Cri, Cre 50 pF 
Feedback network 

Rr 1 GQ Cr ~ 3 pF 

Rk 300 kQ Ck ~ 0.2 pF 

Ce 10 nF 

Signal source circuit 

Ry >1 MQ Ct ~ 0.5 pF 
Note: + indicates the variable resistance range. Without specification, 
the default value after + is 0. 


i.e., Ca © 26 pF in 30 Hz<< f < 3 MHz. The input resistance of Pre-Amp Ra is the 
gate-source resistance of HEMT, and Ra > 10 TQ, so it can be considered as infinity. 

The power supply of Pre-Amp is shown in the dashed box (a2) in Fig/1] A precision 
voltage reference (PVR) LM4050-10 is used in the power supply to provide a constant 
voltage of 10 Volt with typical temperature coefficient of +40 ppm/°C. The noises from 
PVR are eliminated by C1, C2, and Rr. 

For a differential amplifier, there is a feedback between the two brenches of the transis- 
tors, so the hgih common mode rejection ratio (CMRR) can be realized for the differential 
amplifier. For the differential amplifier as Pre-Amp in Ref.|3], if a resistor is instead of 
HEMT (H2) with the same DC resistance, Inv-Amp may becomme a non-inverting am- 
plifier as f — 0. For the design of Pre-Amp in this work, there is no feedback between 
the HEMT brench and Ry brench, therefore Inv-Amp in this work is always an inverting 
amplifier. 


2.2 Design of Post-Amp and composition of Inv-Amp 


In Figfi] the Post-Amp circuit in the dashed box (b) is the same in structure as that in 
Ref. [B]. There is a commercial operational amplifier (OPA) in the circuit, called as Rear- 
OPA, which is the OPA with high gain-bandwidth-product, such as THS4021, OPA657, 
and LMH6624, etc. The Rear-OPA in this work is THS4021 [18]. The Post-Amp parame- 
ters are listed in Table[i] Ra and C, are the input resistance and capacitance of Rear-OPA 
respectively. The feedback resistor R¢ is connected to output O of Rear-OPA and its in- 
verting input. Rr, is 200 kQ, and it is placed in the cryogenic zone. Post-Amp also contains 
two cables that connect Q; and Q» in room temperture zone to outputs O4 and Oz of 
Pre-Amp in cryogenic zone respectively. The capacitance of the two cables is Ci; and Cig 
respectively, and they may vary from 50 pF to 150 pF. In this work, Ci, = Ci2 = Ci = 150 
pF is assumed. 

A 1 meter cable with characteristic impedance of 50 Q has a distributed capacitance of 
100 pF and a distributed inductance of 250 nH, and its distributed LC structure results in 
several resonant gain peaks in [50 MHz, 1 GHz] [19]. Cascade Pre-Amp and Post-Amp to 
form Inv-Amp. In Inv-Amp, Post-Amp and Pre-Amp are connected by the cables (Cable 
OQ, and Cable O2Q,), and the resonant gain peaks may cause the self-oscillations, as 
some unexpected electromagnetic couplings may exist among the outputs and the inputs 
of the amplifiers. Two low-pass filters (Re = Reg = Re = 100 Q and Cy = Cr2 = Cr = 50 
pF) with the upper cut-off frequency of 30 MHz are added at the Rear-OPA inputs, these 
resonant gain peaks can be greatly reduced, thus avoiding the self oscillations. 

For the AC signal, the voltage gain of Inv-Amp is aa = aa(f). aa(f) can be expressed 
as 

aA = Ayp Ayr. (2.3) 


By the nodal analysis method, aa can be obtained. And then, with Ayp expressed by 


Eq. (2.1), Ayr can be obtained by Eq. (2.3). In max{gm/(27Cg),1/(27RsCs)} << f <3 
MHz (i.e. 30 Hz << f < 3 MHz with the parameters listed in Table [1}, 


Re 1 + j2r f RyLCir , 1 
Rup 1+j2rfRaCr 1+ +j2rf iren j 


aa Ry, 


Avr xX (2.4) 


where Ry, = Ry RL/(Ry +R), Cir = Ci+Cr, and aa is the voltage gain of the Rear-OPA. 
In this work, Rẹ >> Ry, and Ra >> RyL. And, a, can be approximately expressed as 
da = aa0o/(1 + j f/f») in (0, 40 MHz] for THS4021. 

With the parameters in Table |1| the Inv-Amp performances are simulated by TINA- 
TI [20]. Fig] show the curves of |aa(f)|ap and Z (aa(f)) simulated by TINA-TI. By the 
nodal analysis method, aa(f) is calculated by the equations established with all compo- 
nents listed in Table |1| [21], which is also shown in Fig[2} The curves of jaa(f)|ag and 
Z (aa(f)) calculated by the nodal analysis method are identical with the simulated ones 
respectively. |aa(f)|ag and Z(aa(f)) are also calculated by Eqs. (2.1), (2.3), and (2.4), 
and the calculated curves in Fig [2] are almost identical with the simulated ones respectively 
in [3 kHz, 10 MHz], which verifies the correctness of Eqs. (2.1), (2.3), and (2.4). 
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Figure 2: The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|ap, 
and the dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI 
simulation results, the green curves are the calculated results with the nodal analysis 
method, and the black curves are the calculated results obtained by Eqs. (2.1), (2.3), and 


EJ. 


2.3 Frequency compensation of feedback loop 


In order to increase the bandwidth of the CryoSTM-TIA, for the high feedback resistor 
Rp with parasitic capacitance Cp, frequency compensation must be used in the feedback 
loop [I| [22]. In Figi] the compensated feedback network is shown in the dashed box 
(c). Taking Ce equalto kCp, where k is above 10°, adjust Ry equalto Rp/k, realizing 
R&C. = RpCp. The output voltage of TIA as V, generates the current Íp flowing to input 
G of TIA, so 

Vo JO Ry + Rf z Rp 

İp 1+j2mfRyCy 1+ j2rfRyCy’ 


where Cy is the parasitic capacitance of Ry, 22]. Zp(f) can be considered as the 
impedance of the feedback network. In Ref.[22], it has been achieved in experiments to 
broaden the bandwidth of the feedback network with the very high feedback resistor Rp 
of 10 GQ to MHz. In (0, 1 MHz], with the parameters listed in Table |1| |Zp(f)| = 
Rp/|1 + j2r f RkCr| > Rp/1.008 and |Zp(f)| < Rp, so it can be considered that Zp(f) is 
equal to Rp. 

In order to verify the above frequency compensation method, the following experiment 
is done. A resistor with 500 MQ(+25 ppm/°C) at room temperature is chosen as Rp, which 
is 1.16 GO at 4.2 K. Rx is a 390 kQ resistor in series with a0 ~ 20 kQ potentiometer, and Ce 
is a COG ceramic capacitor of 10 nF. In order to easily adjust the frequency compensation, 
a 3 pF capacitor is paralleled with the Rp, so that the Cp is about 3.36 pF. Fig.3 shows 
the experimental results for the frequency compensation, in which |Zp| = Rp/V2 at 905 
kHz, and (Zp) as the phase of Zp is —31.5° at 905 kHz. 


Zp(f) = 


2.4 Circuit stability of the proposed CryoSTM-TIA 


In Fig] the signal source circuit is shown in the dashed box (d), and its parameters are 
shown in Table The differential resistance of the tip-sample tunnel junction (TJ) is 
Ry, which is limited to no less than 107? Rp in this design. The capacitance of TJ is Cy, 
which is in parallel with Ry. And, Cy is estimated as several fF [3]. C = Ca + C1 + Cy in 


p= [Zr Rrlas 
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Figure 3: The experimental results for the frequency compensation. At 4.2 K, Rp is 1.16 
GQ and Cf is estimated as 3.36 pF. As R, = 380 kQ and Ce = 10 nF, the black curve 
shows that the bandwidth of |Zp/Rplap is broadened to 905 kHz, and the red curve shows 
Z(Zp) as the phase of Zp is —31.5° at 905 kHz. 


this work. Cy is at least two orders of magnitude less than Ca + C1, so it can be ignored 
in C and C ~% Ca + Cy. The DC bias & modulated signal voltage source is denoted as 
BMS, which provides the DC bias V; and AC sinusoidal modulated signal voltage V; for the 
CryoSTM-TIA. In the following simulation, C{ is always taken as 0.5 pF. TIA connects the 
signal source circuit to form the CryoSTM-TIA. According to the parameters in Table [I] 
the performances of the CryoSTM-TIA can be simulated by TINA-TI. 

The loop gain Tj, of the proposed CryoSTM-TIA is 


Ti(f) = aa (PEE) = aa (PU/C, 


in which 6(f) is the feedback factor, and its reciprocal is, 
1/B(f) = 1+ Zp[l/Ry + 1/Ra + j2af(Ca + C1)]. (2.5) 


By Eq. (2.5), |1/8(f)| and Z (1/6(f)) can be calculated. Fig/4{a) shows the calculated 
results of |1/6(f)|ag and < (1/8(f)) with Ry = +00, and Fig/4{b) shows those with 
Ry = 1 MQ. |aalap and Z (aq) are also shown in Fig/4] Both figures show |TL]aB = 
axles = [1/90 as < <10-0B in. f > 465 ke ond 7) = 2 (aa) 7 1/8) = 
Z (aq) — 180° — Z (1/8) > —130° in f < 465 kHz. Therefore, the CryoSTM-TIA is stable 
with gain margin more than 10 dB and phase margin more than 50°. 


2.5 Voltage gain and transimpedance gain of the proposed CryoSTM- 
TIA 


With the compensated feedback network mentioned in Sect 2.3} it can be considered that 
Zp is equal to Rp in (0, 1 MHz]. Considering the TJ capacitance Cy, the TJ impedance 
should be Zj = Ry/(1+ j2afRjCy). And, C ~ Ca + Cj. As the AC input voltage V; is 
applied by BMS, the output voltage of the CryoSTM-TIA is Vo, and the voltage gain of 
the CryoSTM-TIA is Ay = Vo/V;. In (0, 1 MHz], by the nodal analysis method, Ay is 


Rp 1 
Ay © . : (2.6) 
x 1 R R ; Rp(Ca+C}) 
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Figure 4: TINA-TI simulation results for the voltage gain Inv-Amp a4, and 1/8(f) (a) with 
Ry = +00, (b) with Ry = 1 MQ. Both figures show |Tilap = |aalap — |1/8(f)lap < —10 
dB in f > 465 kHz and Z (TŁ) = Z(—aa) — < (1/8) = Z (aa) — 180° — Z (1/8) > —130° 
in f < 465 kHz. Hence, the CryoSTM-TIA is stable enough. 


Setting V; = 0 and applying a sinusoidal current source Å in parallel with TJ, the 
output voltage V, is generated at the output of the CryoSTM-TIA. A; = Vo Ji is called 
as the transimpedance gain of the CryoSTM-TIA. In (0, 1 MHz], A; is 


Rp 
Ai (2.7) 
1 R R R C wen 
1 aa TAR TAN jenf a i 


Disconnecting TIA with the signal source circuit, and applying a sinusoidal current 
source [yp into the input of TIA, the output voltage Vor is generated at the output of 
TIA. Ait = Vor/ Tyr is called as the transimpedance gain of TIA. In (0, 1 MHz], Air is 


Considering Ra >> Rp and |aa| >> 1, Air is 
Rr 


Ait xX 
TN TL jf Ea’ 


(2.8) 
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Figure 5: TINA-TI simulation results for Air/Rp, Air is the transimpedance gain of TIA 
in the proposed CryoSTM-TIA. |Air(fnr)/Rrlag = —3 dB and Z (Ait (fnr)/Rr) = 116° 
at fot = 201 kHz. 


In Figi] |Air/Relap and Z (Air/Rp) calculated by Eq. (2.8) with the parameters in 
Table |1| is consistent with their TINA-TI simulation results in [10 mHz, 1 MHz], which 
verifies the correctness of Eq. (2.8) for Air. The upper cut-off frequency of Aip(f)/Rp, ie. 
its —3 dB frequency is far = 201 kHz. Comparing Eq. (2.7) with Eq. (2.8), as Ca +C]+Cy & 
Ca +Cr © Ca and Rj > 107?Rp, the upper cut-off frequency of the CryoSTM-TIA is 
approximately equal to fpr. 


2.6 Transient response of the proposed CryoSTM-TIA 


The simulation results for the transient response of the proposed CryoSTM-TIA are shown 
in Supplemental file 3 [24]. For the CryoSTM-TIA, the time taken from adding the input 
step signal voltage to the output response stably within a certain error is called transient 
response time tr. For the simulations, a resistor with a constant resistance Ro is instead 
of the tip-sample junction. When Ro > 30 MQ, for 0.1% error, tp < 5 us. When 5 
MQ < Ro < 30 MQ, for 1% error, tr < 5 us. When 1 MQ < Ro < 5 MQ, for 3% error, 
t < 5 ps. 


3 Inherent noise of the proposed CryoSTM-TIA 


For the circuit of the proposed CryoSTM-TIA shown in Fig[i] the differential equivalent 
circuit with all noise sources is used to calculate its equivalent input noise. The details for 
the noise calculations are shown in Supplemental file 4 [25]. 


3.1 Equivalent input voltage noise and equivalent input current noise of 
Inv-Amp 


The equivalent input noise voltage and equivalent input noise current of HEMT are de- 
noted as ey and iy respectively. The resistors R1, R2, and Rẹ are in the cryogenic zone of 


4.2 K, and their noises in f > 1 kHz are thermal noise, which can be neglected [3]. Rp is 
in the cryogenic zone of 4.2 K, the thermal noise voltage of the resistor Ry is ery. The 
thermal noise voltage of the resistor R.1 is e1, and that of Re2 is eg. The equivalent input 
noise voltage and equivalent input noise current of the Rear-OPA are denoted as ea and 
ia respectively. These noise sources are independent. The equivalent input noise voltage 
and equivalent input noise current of Inv-Amp are denoted as ea and i, respectively. e3 


is the equivalent input noise voltage PSD of Inv-Amp, 74 is its equivalent input noise 
current PSD, eai% is its equivalent input noise voltage-current PSD, and i,e% is its equiv- 
alent input noise current-voltage PSD. By the nodal analysis method and Wiener-Sinchin 
theorem, ignoring the minor terms, it is obtained that [25] 


Z= Riu, ehu | arate | (1 | =) ia (3.1) 
R? Ae. a Ra] Gh 
a= if + + (27 pp (Fug, +a+9+4) 
Ap \ R? 
R 22 
2 HL a 
+27 f) (Con + Ca + HEC) 9, (3.2) 
— Rin, 
cath R = -ir SE (FMR ++ G+) 
vP H 
: RyL Ry \ i2 
— j2 Catat =Q {es a 3.3 
jarf (Cpr + Ca + 2 a) (1+ SE) B (3.3) 


According to the parameters listed in Table |1| for Inv-Amp in this work, e = 0.28 
(nV)?/Hz and i% = 0.14 (fA)?/Hz at f = 10 kHz, and ei = = 0.09 (nV)? /Hz and a = 4.9 
(fA)?/Hz at f = 100 kHz. For the Macro-OPA in Ref.[3], J, e4 = = 0.5 (aV)? /Hz and ¿å = 0.8 
(fA)?/Hz at f = 10 kHz, and e = = 0.14 (nV)?/Hz and i = = 20 (fA)? /Hz at f = 100 kHz. 
The resistor Ry in Inv-Amp in this work replaces the transistor H2 in the Macro-OPA in 

Ref. [B]. The noise generated by Ry is 2 orders of magnitude smaller than that generated 
by H2 in the Macro-OPA in Ref.[8], i.e. ey << Apez. Therefore, the equivalent input 
noises of Inv-Amp in this work are much lower than those of the Macro-OPA in Ref. [3]. 

In Eqs. (8-1), (3.2), and (3.3), e2 gR / RẸ is one order of magnitude smaller than e. 
In Eq. (8-1), (1 + Rut /Ra)° i2/92, is 2 orders of magnitude smaller than e2,. In Eq. (3.2), 
(Ces + Cga + Ca RuL/ Ra)? i2/92, is one order of magnitude smaller than C?e2/A2p. In 
Eq. (3.3), (Ces + Cga + Ca Rut /Ra) (1 + Rut /Ra) i2/ g2 is one order of magnitude smaller 


than Cye2/A?,. Further ignoring the minor terms in Eqs.(3.1), (3.2), and (3.3), 
a vP 


eh © eh + (2 + ef + 63)/Atp, (3.4) 
TR m i + (2m f)?OX (€2 + ef + e2) /Ap, (3.5) 
eath © (aeh )* = —j2nfCa(e2 + e? + €3)/Arp. (3.6) 
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3.2 Equivalent input voltage noise and equivalent input current noise of 
the proposed TIA 


The equivalent input noise voltage PSD of TIA is denoted as e2, its equivalent input 


noise current PSD is i2, its equivalent input noise voltage-current PSD is epi, and its 
equivalent input noise current-voltage PSD is ireț- They are 


ep =e ~ ef + (e2 +e} + €5)/Adp, (3.7) 


i} = 4kpT/Rp + eh /Rp +i 
etare 
Ap RE 
+42, + (2m f)?CR (e2 + ef + e3)/A?p, (3.8) 


~ 4kpT/Rp + €2,/R2 + 


eri} = (ire )* = e /Rp + eai 


-j2nfCa (e? + e? + e3)/A?p. (3.9) 
3.3 Equivalent input current noise of the proposed CryoSTM-TIA 
The equivalent input noise current PSD of the proposed CryoSTM-TIA is 


2 _— + [1/R5 + (27 f) 202] e2 e, 
+(1/Ry a. (1/Ry eee (3.10) 


where Cj = C1 + Cy. Putting Eqs.(3.7| , and (3.9) into Eq.(3.10} (3.10), 


iZ =i, + 4kpT/Rp + (1/Ry +1/Rp)? = +(2@+et+e) /A2>| 


+ (Inf)? lc? Tyee + C7(e2 + ef + )/Azp| , (3.11) 


i2 =i + 4kgT/Re + (1/Ry + 1/Rp) & + By, (3.12) 


where iĝ = (27 f)? eae + 2C',C13(e2 + e? + €2)/A2, 

For the proposed CryoSTM-TIA, Rp = 1 GO, Ca = 26 pF, Ci = 0.5 Be and Cy = 10 
fF. Rp and TJ are in the cryogenic zone at 4.2 K. As Ry = 1 MQ, i2 and its four 
components are listed in Table [2] The noise components of the CryoSTM-TIA in Ref. [3] 
are also shown in Table The CryoSTM-TIA proposed in this work has an equivalent 
input noise current PSD of 0.62 (fA)? /Hz at 10 kHz and 5.3 (£A)? /Hz at 100 kHz, which 
is much lower than that in Ref.[8]. 

For the Epi-wafer of the CNRS-HEMT, as Si is doped in Al,Ga,_,As, each dopant 
gives rise to two types of electronic states: a shallow and delocalized donor level associ- 
ated with the normal substitutional site configuration, called as shallow donor state; a 
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Table 2: Noise components of CryoSTM-TIAs 


Pre-Amp type Single HEMT Amp Pre-Amp in Ref. [3] 
f (kHz) 10 100 10 100 
e? ((nV)*/Hz) 0.28 0.09 0.5 0.14 

Unit for the following terms is (fA)?/Hz 

in 0.14 4.9 0.8 20 
4kpT/ Rr 0.2 0.2 0.2 0.2 
B 0.0033 0.15 0.03 0.7 
(ae tei 0.28 0.09 0.5 0.14 
i2, as total 0.62 5.3 1.5 21 


in 


more localized acceptor level, known as DX~ center, arising from a lattice distortion. In 
Al,Ga ,_,As, as x > 0.22, DXT as a bound state is more stable than the shallow donor 
state. And, x = 0.37 for CNRS-HEMT in this work. There is a large repulsion barrier 
Ecap in the transition from the shallow donor to DX~, so DX- is frozen as T < 120 K 
[27]. According to the work on the noise mechanism of the CNRS-HEMT [5], “bias- 
cooling method” can be used to control the density of the frozen DX- centers in the 
doping area in the HEMT. This unique method, by which the in-situ modification of the 
HEMT structure is realized, can be used to reduce the noise generated by the HEMT. By 
cooling the HEMT with the positive gate-source voltage, more electrons in the HEMT are 
frozen into the DX” centers. At low temperatures, the absolute value of the gate-source 
voltage required to maintain the HEMT at the ideal operating point is reduced, so that 
the gate leakage current is reduced, and thus the low frequency noises generated by the 
gate leakage current is reduced [5]. By the means, the inherent noise of the CryoSTM-TIA 
can be reduced. Moreover, if the HEMT used for the CryoSTM-TIA is not CNRS-HEMT, 
or its noise performances are not as good as those listed in Table [I] it can be remedied 
by “bias-cooling method”. This approach is very convenient for cases with only a single 
HEMT in Pre-Amp. 


4 Operating state adjustment and DC tunneling current 
measurements 


For the CryoSTM-TIA in Fig[i] along with the circuit parameters listed in Table |1} its 
operating state can be adjusted as follows. (1) Disconnect Pre-Amp from Post-Amp and 
ground the input of Pre-Amp (i.e. the Gate G of HEMT). Keep Rp to 0. Adjust Rs, Rp1, 
and Rt, to achieve HEMT at the ideal operating point (Vag = 100 mV, Ig; = 1 mA) and 
the potentials at O} and Og equal. That is to say, Rg; and Ry, are adjusted to realize 
Ry = Rg + Vas/Ias = Rs +100 Q. (2) Cascade Pre-Amp and Post-Amp to form Inv-Amp, 
and input G is still grounded. Adjust Rp and Ri; to realize the DC output voltage Vom 
of Inv-Amp to 0 and HEMT still at the ideal operating point. (3) Connect output O and 
input G of Inv-Amp with the feedback network and disconnect input G from ground, to 
form TIA. Input G is still not connected with the signal source circuit, so the potential 
of input G and output O of Inv-Amp is still 0, since the input resistance of the transistor 
HEMT Ra can be considered as infinity. In the feedback network, Rp + Rk + Rp = 1 
GQ. (4) Connect the signal source circuit to TIA, to form the CryoSTM-TIA. 
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When f — 0, aa > aao, and aao is the DC voltage gain of Inv-Amp. |aaol|ap ~ 
(gmRe/gaRu)ap ~ 94.5 dB, which is consistent with the simulation result |aao|azn = 92 
dB shown in Fig P| As the DC bias V; is applied by BMS, the DC resistance of TJ is 
R, and the potential at input G is Vg, and the output voltage of the CryoSTM-TIA is 
Vo. Obviously, aagVe = Vo, and (Vi — Va)/R = (Ve — Vo)/Rr. The DC bias on TJ is 
V =V,— Vg, and the DC tunneling current J = (V; — Va )/R. I; as an approximate value 
of I is 

Is = —Vo/Rr. (4.1) 


And, the relative error is obtained as 
Er = |I; — I|/|I| = 1/(1 — aao). 


As |aao|ap © 92 dB, Er < 30 ppm, consistent with the simulation results [28]. Since 
aaoVa = Vo, the DC bias on TJ is 


VV a. (4.2) 


By Eq.(£1), Vo = —Rrls ~ —RpI = —RpV/R. By Eq.(42), Vi ~ V — ReV/(Raao). As 
the minimum of Ry is not less than 1 MQ, R > 1 MQ, so |Rp/(Raao)| < 1/40. Thus, 
V x Vi. With the measured V, I ~ J, and V & V; can be obtained. Therefore, the 
scanning tunneling current spectra I = I(V) (V €[Vi, Vu) can be obtained, where Vg 
and Vy are the lower limit and upper limit voltages for the measurements respectively. 

The drift of Vom can be considered as the amplification of the drift of the Inv-Amp 
input offset voltage Vos, where Vom = GagVos. For TIA, the output drift voltage is 
approxmate to Vog. The estimated drift of |Vog| is less than 17 wV/°C [29]. The power of 
the power supply of Pre-Amp is less than 40 mW. The temperature control system based 
on the TEC devices can be used to control the temperature fluctuations of the power 
supply and those of Rear-OPA within 0.01 °C, so the fluctuations of Vog within 170 nV 
is guaranteed, i.e. the TIA output fluctuations within 170 nV. 

In order to verify the correctness of the above design methods for TIA, we do the fol- 
lowing work shown in Ref.|31]. The N-Channel JFET (SST4393-T1) is chosen to replace 
CNRS-HEMT, and TIA is designed and fabricated with the same circuit design method 
and the same circuit topology. SST4393-T1 can work at 77 K. We measured the electrical 
performances and noise performances of the JFET-based TIA at 77 K [BI]. The experi- 
mental results are basically consistent with the calculated results and simulated results, 
which shows that the design methods for TIA in this work is correct. 


5 Applications for the proposed CryoSTM-TIA in spectra 
measurements 

For most applications, the modulus of the transimpedance gain of the CryoSTM-TIA 
|Ai(f)| should be measured firstly. Aj(f) is expressed as Eq.(2.7). Since |aa(f)|ap > 90 
dB in [1 kHz, 300 kHz] as shown in Fig|2} when Ry > 107° Rp, |Ai(f)| in [1 kHz, 300 kHz] 
can be approximately expressed as 

|Ai(f)| = Re/|1 — j2r f Rr (Ca + C1)/aa (f). (5.1) 
How to measure it has been described in Ref. [3]. 
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5.1 Measurements of scanning tunneling differential conductance spec- 
tra by the proposed CryoSTM-TIA 


The differential conductance of TJ Gy = 1/Ry is the function of the voltage V applied 
to TJ. As the frequency f of the modulated signal voltage Vi provided by BMS is low 
enough, such as f < 1 kHz, by Eq.(2.6), 

Rp 1 


Ay & . 
RV) 1- 


Rr i 
aa (f)R3(V) 


Therefore, Ry(V) ~ |-1/Ay + 1/aa(f)] Rr. When the measured |A,| < 1000, R3(V) ~ 
Rp /|Ay|, since |aa(f)|ag > 85 dB in f < 300 kHz as shown in Fig [2] With the measured 
|Ay| = |Vo|/|Vi|, the differential conductance spectra G3(V) = 1/R;(V) (V €[Vz, Va]) can 
be obtained. 

Increasing the frequency f of the modulated signal can speed up the scanning tunneling 
differential conductance spectra measurements. In [1 kHz, 300 kHz], as Ry > 1 MQ, by 
Eq. (2.6), 

(AIS po fir 

. IZa(f)| |1— g2afRr(Ca + Cr)/aa(f)| 
Therefore, 1/|Z3(f)| ~ |Av(f)|/|Ai(f)|, where | Ai(f)| is shown by Eq. (5-1). And, |Ay(f)| 
and |A;(f)| can be measured. 1/|Z3(f)| = 1/R3 + (2nfCy)* can be obtained. Selecting 
two different frequencies fı and fz in [1 kHz, 300 kHz], |Z3(f1)| and |Z3(f2)| are obtained. 
Ry and Cy can be solved out from the measured |Z 3(f1)| and |Z3(f2)|. 

Since the inherent noise of the CryoSTM-TIA is very small, the amplitude of the 
modulated signal voltage V, can be very small, so that the energy resolution for STS is 
much improved. 


5.2 Measurements of scanning tunneling shot noise spectra by the pro- 
posed CryoSTM-TIA 


The measurement method for the tunneling shot noise spectra by the proposed CryoSTM- 
TIA is basically the same as that introduced by Ref.[3], and it is only briefly described 
here. 

Before approaching the tip to the sample in the CryoSTM, Rj can be considered as 
infinity and Cy as 0. In this case, the output noise voltage PSD of the CryoSTM-TIA 
Ssu(f) can be measured, and the equivalent input noise current PSD i?(f) is 


P(f) = Sul PAA, (5.2) 
where |Aj;(f)| is obtained by Eq. (5.1). 


To measure the scanning tunneling shot noise spectra (STSNS) of a quantum system by 
the CryoSTM-TIA, the distance between the tip and sample is adjusted, and the interval 
(VL, Vu] is selected, so that the shot noise measurements at V € Dy = {V| V% < V < 
Vu,G3(V) < 1 uS} are needed to study the physical properties of the quantum system 
[3]. In Dy, the output noise voltage PSD Ssum(f, V) can be measured. The tunneling 
noise current PSD is denoted as S}(f,V). The equivalent input noise current PSD of the 
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CryoSTM-TIA is the function of Ry, so it is the function of V and denoted as i2 ( f,V). 
For Ssum( f, vV), 


Ssum(f,V) = [IEN EEN] AAP. (5.3) 
With Eq. (5.2) and (5.3), it is obtained that 


Ssum( f, V) _ Ssul f) 
|Ai(f)|? 


where 6(f,V) = 2,(f,V) — #2(f). From Eq. (3.12), 


2 a , 4kpT | eA 2 | 7200 ateta 


? "In 


Iô(f, V)| is smaller than 0.3 (fA)?/Hz in [10 kHz, 100 kHz) and 0.1 (£A)? /Hz in[100 kHz, 
200 kHz]. ô(f,V) can be neglected compared with S1(f, V), so long as the minimum of 
G3(V) (V € Dy) is not too small. Therefore, S1(f, V) can be obtained as 


Ssum ald — Ssu 
SIV) m 


by the measured Ssum(f, V), Ssu(f), and |Ai(f)|. And then, the SNTNS Sj,(J) = 2Fe|J| 
can be extracted from S{(f,V). Here, J = I(V) (V € Dy), and I is the DC tunneling 
current as the bias V is applied. 

Compared with the CryoSTM-TIA in Ref.[3], the CryoSTM-TIA proposed in this work 
has the same transimpedance gain, the bandwidth is above 200 kHz, and its inherent noise 
is much lower. The inherent noise of the CryoTM-TIA in this work is only 5.3 (fA)?/Hz 
at 100 kHz and its 6(f,V) < 0.1 (fA)?/Hz at 100 kHz, while that of the apparatus in 
Ref.[3] is 21 (fA)?/Hz at 100 kHz and 6(f,V) < 0.28 (fA)?/Hz at 100 kHz. Therefore, 
as investigating novel quantum states of various quantum systems, the measurements 
performed with the apparatus proposed in this work is more accurate. 

For example, as investigating the existence of MBS in a magnetic flux vortex of the 
iron superconductor in CryoSTM [4], the tunnel junction resistance must be large enough, 
since the occurrence of incoherent Andreev reflection can be ruled out by weak tunnel 
coupling conditions [32]. In Ref.[4], the tunneling current I is quite low (see Fig.3(a) in 
Ref.{4] and Fig.s2(a) and (b) in its supplemental file). As the tunnel junction bias V is 
0.3 mV, I is only tens of pA, and the corresponding shot noise may only be only a few 
(fA)?/Hz. For the shot noise measurements in the system in Ref.[4], it is obvious that the 
CryoSTM-TIA proposed in this work is much more effective and accurate than that in 
Ref. 3]. 


And, i? (f, V) is composed by the four parts as shown in Eq. (3.12). According to Table |2| 


(5.5) 


6 Conclusion 
In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning tunneling 


microscope (CryoSTM) is presented. TIA connected with the tip-sample component in 
CryoSTM is called as CryoSTM-TIA. The CryoSTM-TIA in this work has transimpedance 
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gain of 1 GO and bandwidth more than 200 kHz, and its equivalent input noise current 
PSD is only 5.3 (fA)?/Hz at 100 kHz. In this CryoSTM-TIA, a single CNRS-HEMT 
Pre-Amp is instead of the differential Pre-Amp consisting of a pair of CNRS-HEMTs in 
Ref.[B]. The difficulty of matching the identical HEMTs is avoided. And, due to the 
reduction of a noisy HEMT, the apparatus inherent noise is only 1/4 of that in Ref.[3}. 
Furthermore, the inherent noise of the single HEMT in the circuit can be in-situ reduced 
by “bias-cooling method”. With this apparatus, the fast high-energy-resolution scanning 
tunneling spectra measurements can be performed and the very low tunneling shot noise 
of quantum systems can be measured at the atomic scale. This apparatus can be applied 
to investigate novel physical properties of various quantum systems, such as detecting the 
existence of Majorana bound states in the topological quantum systems. 


Acknowledgment 


This work is supported by Open Research Fund Program of the State Key Laboratory of 
Low-Dimensional Quantum Physics (Grant No. KF202212). I acknowledge stimulating 
discussions with Prof. Fang-Hao Liang of Math School in Shandong Univ and Prof. Xing 
Liang of Math School in USTC. 


Declaration of competing interest 


The author declares that they have no known competing financial interests or personal 
relationships that could have appeared to influence the work reported in this paper. 


References 


[1] C.J. Chen, Introduction to scanning tunneling microscopy, Oxford Univ. Press, 
(1993). 


[2] J.F. Ge, M. Ovadia, and J.E. Hoffman, Achieving low noise in scan- 
ning tunneling spectroscopy, Rev. Sci. Instrum. 90 (2019) 101401, 


https://doi.org/10.1063/1.5111989 


[3] Y.X. Liang, Low-noise large-bandwidth tansimpedance amplifier for measuring scan- 
ning tunneling shot noise spectra in cryogenic STM and its applications, Ultrami- 


croscopy 234 (2022) 13466, https: //doi.org/10.1016/j.ultramic.2022.113466 


[4] D.F. Wang, L.Y. Kong, P. Fan, H. Chen, S.Y. Zhu, W.Y. Liu, L. Cao, Y.J. Sun, S.X. 
Du, J. Schneeloch, R.D. Zhong, G.D. Gu, L. Fu, H. Ding, and H.J. Gao, Evidence 
for Majorana bound states in an iron-based superconductor, Science 362 (2018) 333, 


https: //doi.org/10.1126/science.aaol797 


[5] Y.X. Liang, Q. Dong, M.C. Cheng, U. Gennser, A. Cavanna, and Y. Jin, In- 
sight into low frequency noise induced by gate leakage current in AlGaAsGaAs 
high electron mobility transistors at 4.2 K, Appl. Phys. Lett. 99 (2011) 113505, 


https: //doi.org/10.1063/1.3637054 


[6] 


[8] 


[9] 


10 


11 


12 


[15] 


[16] 


Y. Jin, Q. Dong, A. Cavanna, U. Gennser, L. Couraud, and C. Ulysse, Ultra-low noise 
HEMTs for deep cryogenic low-frequency and highimpedance readout electronics, 
12th IEEE International Conference on Solid-State and Integrated Circuit Technology 
(ICSICT) (2014). 


K.M. Bastiaans, D. Cho, T. Benschop, I. Battisti, Y. Huang, M.S. Golden, 
Q. Dong, Y. Jin, J. Zaanen, and M.P. Allan, Charge trapping and super- 
Poissonian noise centres in a cuprate superconductor, Nat. Phys. 14 (2018) 1183, 


https: //doi.org/10.1038/s51567-018-0300-z 


F. Massee, Y.K. Huang, M.S. Golden, and M. Aprili, Noisy defects in 
the high-T, superconductor BigSrgCaCu,O0g;2, Nat. Commun. 10 (2019) 544, 


https: / /doi.org/10.1038/s51467-019-08518-1 


C.J. Bolech and E. Demler, Observing majorana bound states in p-wave superconduc- 
tors: using noise measurements in tunneling experiments, Phys. Rev. Lett. 98 (2007) 


237002, https: //doi.org/10.1103/PhysRevLett.98.237002 


A. Golub and B. Horovitz, Shot noise in a Majorana fermion chain, Phys. Rev. B 83 


(2011)153415, |https://doi.org/10.1103/PhysRevB.83.153415 


H. Soller and A. Komnik, Charge transfer statistics of transport through Majorana 


bound states, Physica E 63 (2014) 99, https://doi.org/10.1016/j.physe.2014.05.020 


C.W.J. Beenakker and D.O. Oriekhov, Shot noise distinguishes Majorana fermions 
from vortices injected in the edge mode of a chiral p-wave superconductor, SciPost 
Phys. 9 (2020), 080, https://doi.org/ 10.21468/SciPostPhys.9.5.080 


K.M. Bastiaans, D. Cho, D. Chatzopoulos, M. Leeuwenhoek, C. Koks, and M.P. 
Allan, Imaging doubled shot noise in a Josephson scanning tunneling microscope, 


Phys. Rev. B 100 (2019) 104506, https: //doi.org/10.1103/PhysRevB.100.104506 


K.M. Bastiaans, D. Chatzopoulos, J.F. Ge, D. Cho, W.O. Tromp, J.M. van 
Ruitenbeek, M.H. Fischer, P.J. de Visser, D.J. Thoen, E.F.C. Driessen, T.M. 
Klapwijk, and M.P. Allan, Direct evidence for Cooper pairing without a spec- 
tral gap in a disordered superconductor above Te, Science 374 (2021) 608. 


https: //doi.org/10.1126/science.abe3987 


S. Cocklin and D.K. Morr, Scanning tunneling  shot-noise spec- 
troscopy in Kondo systems, Phys. Rev. B 100 (2019) 125146, 


https: //doi.org/10.1103/PhysRevB.100.125146 


Webpage of HF-STM-broadband cryogenic buffer amplifier made by Stahl-electronics, 
https: //www.stahl-electronics.com/products/amp/hfstm/, In its data sheet (pdf), 
the input capacitance of this amplifier, including the input capacitance of the HEMT 
and the capacitance of the wire between the tip and the HEMT input, is only 2.7~3.1 
pF. The input capacitance of the HEMT C% is 2.5 pF at least [6]. Therefore, Cy < 0.5 
pF is reasonable. 


17 


17 


18 
19 


20 


21 
22 


23 


24 


25 


26 


[27] 


28 
29 


30 


31 


[82] 


Data sheet of LM4050-10, https://www.ti.com.cn/document-viewer /cn/lm4050- 
n/datasheet 


Webpage of THS4021 OPA, https://www.ti.com/product /THS4021 


Supplemental file 1. 


TINA-TI is SPICE-based analog simulation program produced by Texas Instruments 


Inc., https://www.ti.com/tool/TINA-TI 


Supplemental file 2. 
B. Michel, L. Novotny, and U. Dürig, Ultramicroscopy 42-44 (1992) 1647, 


S. Franco, Design with Operational Amplifiers and Analog Integrated 
Circuits,McGraw-Hill Companies, Inc., (2002). 


Supplemental file 3. 
Supplemental file 4. 


D.J. Chadi and K.J Chang, Energetics of  DX-center formation 
in GaAs and AlxGal-xAs alloys, Phys. Rev. B 39 (1989) 10063, 
https://doi.org/10.1103/PhysRevB.39.10063 


W.F. Li, Y.X. Liang, Y. Jin, and J.H. Wei, Bipolaron mechanism 
of DX center in Al,Ga,_,As: Si, Acta. Phys. Sin. 59 (2010) 8850, 
https: //doi.org/10.7498/aps.59.8850 


Supplemental file 5. 


Supplemental file 6. 


A suitable TEC device, |https://datasheets.maximintegrated.com/en/ds/MAX1978- 


MAX1979.pd 


Y.X. Liang, Low-noise  large-bandwidth high-gain transimpedance am- 
plifier for cryogenic STM at 77 K, submitted to  Ultramicrscopy, 
https: //chinaxiv.org/abs/202404.00274 


L.Y. Kong and H. Ding, Emergent vortex Majorana zero mode 
in iron-based superconductors, Acta. Phys. Sin., 69 (2020) 110301, 
https: //doi.org/10.7498 /aps.69.20200717 


Supplemental file 1: Voltage gain peaks for cables 
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Figure s1-1 The circuit scheme for measuring the voltage gain peaks for cables. 


A 1 meter cable with characteristic impedance of 50 Q has a distributed capacitance 
of 100 pF and a distributed inductance of 250 nH, which is shown as Fig. s1-1. In Fig. s1- 
1, R=1.5 kQ, the gain is Vo/Vin. The simulation results show its distributed LC structure 
results in several resonant gain peaks in [80 MHz, 1 GHz] shown as Fig. sl-2. As the 
input signal Vin is the thermal noise of R (1.5 kQ) at 300 K, we amplify Vo with 
Amplifier SA-220F5 (gain=200, fr=200 MHz), and then the gain peaks for the cable are 


measured with an oscilloscope. Fig. s1-3 show the measured gain peaks for cable. 
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Figure s1-2 The simulation results of the gain peaks for a 1 meter cable with a distributed capacitance 
of 100 pF and a distributed inductance of 250 nH. 
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Figure s1-3 The input signal Vin is the thermal noise of R (1.5 kQ) at 300 K. The measured results of 
the gain peaks for a 1 meter cable with a distributed capacitance of 100 pF and a distributed 
inductance of 250 nH. 


Supplemental file 2: Voltage gain of Inv-Amp 
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Figure s2-1 Inv-Amp circuit with all components for Eqs.(s2.1-s2.8), and their values are shown in 
Table 1 in Article. 


With the nodal analysis method, Inv-Amp voltage gain aa(f) is calculated by the 
following equations. All components for Eqs.(s2.1-s2.9) are shown in Fig.s2-1, and their values 


are shown in Table 1 in Article. 
(V, -V,)/R, =(V, —V,)/R, +(V, =R +7, /Ry + joC,V, > (s2.1) 


JoCa (P-V) = 8a (V-V.)+ a (P-P) +0, -0 )/R +P P.) Ra FCP, 622 


( -P )/R +% -0 )/R, =V,/R, + j@C P> (82.3) 
(V,-V,)/R, =(V_-%)/R, JE s (s2.4) 
En (VV) + ga (P -8 )=P/R + JOCV, + jC, (V.-V,), (s2.5) 


V,=a,(V,-V), (s2.6) 


(V, ~ V.) /Ra =+joC,V, (s2.7) 


(V,-V.)/R, =+joC,V,, (s2.8) 


a, =V,/V, (s2.9) 
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Figure s2-2 The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|as, and the 
dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI simulation results, the green 
curves are the calculated results with the nodal analysis method (i.e. Eqs.(s2.1-2.9)), and the black 
curves are the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. 


Fig. s2-2 shows the modulus of the Inv-Amp voltage gain |aa(f)|az. The red solid 
curve is the TINA-TI simulation results, and the green solid curve is the calculated results 
with the nodal analysis method. They are identical with each other. The black solid curve 
is the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. In fe[3 kHz, 10 MHz], 


three curves are identical with one another. 


Fig. s2-2 also shows the argument of the Inv-Amp voltage gain Z(aa(f)). The red 
dashed curve is the TINA-TI simulation results, and the green dashed curve is the 
calculated results with the nodal analysis method. They are identical with each other. The 
black dashed curve is the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. In 


Jel kHz, 10 MHz], three curves are consistent with one another. 


Supplemental file 3: Transient response of the proposed CryoSTM-TIA 
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Figure s3-1 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with the 
different Ro for the time interval of 10 us. The dashed curve is the step input signal Vi/Vinst for the 
different Ro. The solid curves are the output response Vo/Vost. Here, Vost and Vinst are their values 
respectively at 300 ms after applying the step input signal Vi, which are shown in Fig.s2-2. Transient 
response time t<5 us. (b) is the zoom-in for the output response Vo/Vost in (a). 
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Figure s3-2 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 


simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with different 


Ro for the time interval of 300 ms. (a) The step input signal Vi/Vinst. (b) The output response Vo/Vost. 


Vo/Vos= 1 in t> 100 ms. (c) is the zoom-in figure of (b). 


A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. Fig.s3-1 shows TINA-TI simulation results for the transient response of the 
CryoSTM-TIA with the different Ro for the time interval of 10 us. The dashed curve is 
the step input signal Vi/Vins, and the solid curves are the output response Vo/Vost for the 
different Ro. Here, Vost and Vinst are their values respectively at 300 ms after applying the 
step input signal Vi. The error is definite as |Vo—Vost |/|Vost|. For the CryoSTM-TIA, the time 
taken from adding the input step signal voltage to the output response stably within a 
certain error is called transient response time f. As shown in Fig.s3-1 and Fig.s3- 
2(b)&(c), with Ro= 1 GQ and 30 MQ, for 0.1% error, t<5 us. With Ro=5 MQ, for 1% 
error, &<5 us. And, With Ro= 1 MQ, for 2.5% error, tr <5 us. 


Supplemental file 4: Noise of the proposed CryoSTM-TIA 


S4.1 Noises of STM-TIA 
For Inv-Amp, the equivalent input noise voltage of Inv-Amp is denoted as ea and its equivalent 


input noise current is ia, and their harmonic components of frequency f are E and Ja respectively. 


E 
i) can be obtained by the nodal analysis method. By Wiener-Khintchine theorem, 


A 


= 
ea Eala 


EE EL 


AA AVA 


can be obtained from l [S4R1, S4R2]. The two matrix elements on 


DEn in 
the main diagonal are the equivalent input noise voltage PSD of Inv-Amp e} and its equivalent 


input noise current PSD i, . The two matrix elements on the sub-diagonal are its equivalent input 


noise voltage-current PSD e,i, and equivalent input noise current-voltage PSD i,e, . 


Inv-Amp is connected to the feedback resistor Rr to form TIA. The equivalent input noise 
voltage of TIA is denoted as er and its equivalent input noise current is ir, and their harmonic 
components of f are Er and Ir respectively. The temperature of the feedback resistance Rr is T, and 


its noise voltage is er and its harmonic component of fis Er. 


(b) 


Figure s4-1 (a) TIA circuit with the input short-circuit containing the equivalent input noise voltage of Inv-Amp 
Inv e4 and its equivalent input noise current ia, and the output noise voltage of eorv ; (b) Noiseless TIA circuit with 
the equivalent input noise voltage of TIA er as the input signal, and the output noise voltage of eorye ; the 


equivalency of the above two circuits means eoty = €otve - 


For TIA, the circuit containing all noise sources with the input short-circuit is shown as 


Fig.s4-1(a), and the output noise is eoty. The noiseless circuit with the equivalent input noise voltage 
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of TIA er as the input signal is shown as Fig.s4-1(b), and the output noise is eotve. For calculating 
the equivalent input noise voltage of TIA, the equations are established on the equivalency of the 
above two circuits, i.€. Coty = @otve. Therefore, by nodal analysis method, 


ET=E,. 


Figure s4-2 (a) TIA circuit with the input open-circuit containing the equivalent input noise current of Inv-Amp 
Inv ea and its equivalent input noise current ia, and the output noise current of eori ; (b) Noiseless TIA circuit with 
the equivalent input noise current of TIA ir as the input signal, and the output noise voltage of eotie ; the 


equivalency of the above two circuits means eoti= €oric - 


For TIA, the circuit containing all noise sources with the input open-circuit is shown as 
Fig.s4-2(a), and the output noise is is eori. The noiseless circuit with the equivalent input noise 
current of TIA ir as the input signal is shown as Fig.s4-2(b), and the output noise is eotie. For 
calculating the equivalent input noise current of TIA, the equations are established on the 
equivalency of the above two circuits, i.€. eoTi = otic. Therefore, by the nodal analysis method, 


Tp=Iy+(Eat EpyRe . 


lun ee 


i : = 
From ee a , by Wiener-Khintchine theorem, a 
LE, hh Le i 


Hence, 


a 
erly 


can be obtained. The two 


matrix elements on the main diagonal are the equivalent input noise voltage PSD of TIA e and its 
equivalent input noise current PSD a The two matrix elements on the sub-diagonal are its 


equivalent input noise voltage-current PSD e,j, and equivalent input noise current-voltage PSD 


* 


iper: 
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= =e (s4.1) 


P =P +4k,T/R,+e2/R2 , (s4.2) 


ei = (ine )* ent, ten /Re (543) 
Here, the noise voltage PSD of Rr is g = 4k, TR; . 
TIA is connected with the signal source circuit to form a STM-TIA. The equivalent input noise 


current PSD of a STM-TIA i can be obtained [S4R1-S4R3] by 
in =i +(YR +(2nf) Cy )er +(I/R, + J2mfCy erie +(/R — j2mfCy ine. (s4.4) 


Here, Cy=Crt+Cj. Putting Eq.(s4.1), ($4.2), and (s4.3) into Eq.(s4.4), ie 


=i? +4k,T/R; +3 (1/R? +1/R? + (2nf)° Ci) 


(84.5) 
ae + j2mfC,)( enix +e R) +(I/R,— j2mMCy )(i e, + eK / Re) 


S4.2 Noises of the proposed CryoSTM-TIA 
S4.2.1 Equivalent input noise voltage and the equivalent input noise current of Inv-Amp 

For Inv-Amp, the equivalent input noise voltage and equivalent input noise current of HEMT 
are denoted as en and in respectively, and their harmonic components of f is are denoted as Ey and 
Iy respectively. gm is the transconductance of HEMT and ga is channel conductance of CH. The 
noise voltage of the resistors Ry is denoted as ern, and the harmonic components of fare denoted as 
Eu. The noise voltage of the resistors R1, R2, and the feedback resistance Rr are denoted as ex1, e12 
and er respectively, and the harmonic components of fare denoted as E11, E12, and Ef respectively. 
The noise voltage of the resistors Rc: and Re2 are denoted as e1, e2 and ef respectively, and the 
harmonic components of f are denoted as E1, E2, and Ef respectively. Here, Ri= Ro= Rv. All noises of 
Ru, Ri, R2, Rei, Re2, and Rs above 1 kHz are thermal noise. Ra= Ri//(1/ga) and Rut = Ru//Ri are 
denoted. The equivalent input noise voltage and equivalent input noise current of Rear-OPA are 
denoted as ea and ia respectively, and the harmonic components of f are denoted as Ea and Ja 
respectively. These noise sources are independent of each other. For Inv-Amp consistent of 
Pre-Amp and Post-Amp, its equivalent input noise voltage and equivalent input noise current as ea 


and ia, the corresponding harmonic components of fare Ea and Ta respectively. 
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Figure s4-3 (a) Inv-Amp equivalent differential circuit with the input short-circuit containing all noise sources, and the 
output noise current of eoay; (b) Noiseless Inv-Amp circuit with the equivalent input noise voltage of Inv-Amp ea as the 
input signal, and the output noise voltage of eoave; the equivalency of the upper two circuits means eoay = oAve. In (a) 


and (b), the triangle OPA is the Rear-OPA. 


The equivalent differential circuit of Inv-Amp containing all noise sources with the input 
short-circuit is shown as Fig.s4-3(a), and the output noise is eoay. The noiseless circuit with the 
equivalent input noise voltage of Inv-Amp ea as the input signal is shown as Fig.s4-3(b), and the 
output noise is eoave. For calculating the equivalent input noise voltage of Inv-Amp, the equations 
are established on eoav = ove. Therefore, by the nodal analysis method, 

Ea ~ Eny + ErnZuv/(Avei Ru) EL (EmRL) + Za(E2/Ri + Ei/Re)/Avei 

+ Ea/Avpj + (E1-E2)/Ave1 ~ (1/8m + Zuu Avrj)ha , (s4.6) 
where Ayp1 = 2m/[1/Ritgaty(2nf/)Ci], Ave) =2m/[1/Ritgaty(2rf\Cir], Zar = Ruv/[1+ j2a)RarCir], and 
Cir =C; +. 

The equivalent differential circuit of Inv-Amp containing all noise sources with the input 
open-circuit is shown as Fig.s4-4(a), and the output noise is eoai. The noiseless circuit with the 


equivalent input noise current of Inv-Amp ia as the input signal is shown as Fig.s4-4(b), and the 
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output noise is eoaie. For calculating the equivalent input noise current of Inv-Amp, the equations are 
established on @oai= odie. Therefore, by the nodal analysis method, 
Ta~Tq t](20f)CarErnZur/(AveiRu) — J(27f/)CsaEui/(gmRi) + j(2nf\CaiZar(Ev2/Ri+Ed/Rs)/Aye 1 

+ j(2nf) (E1-E2)/Avei— j(21f)CajE al Avr j (27) (Csa/gm + CajZu/A vej)la, (s4.7) 
where Cai = Cgst+Cga(1+Ave1) , Csa = CgstCga, and Caj = CystCga(1+Avpj). 


(a) 


Figure s4-4 (a) Inv-Amp equivalent differential circuit with the input open-circuit containing all noise sources, 
and the output noise voltage of eoai; (b) Noiseless Inv-Amp circuit with the equivalent input noise current of the 
Inv-Amp ja as the input signal, and the output noise voltage of @oaie; the equivalency of the upper two circuits 


means €oai= odie. In (a) and (b), the triangle OPA is Rear-OPA. 


For the measured frequency of fm < 300 kHz, 2xfC;<0.3 mS with C=150 pF and 2zfC,<0.1 mS 
with C.=50 pF while 1/Ra = 1/Ri+ga = 2 mS and 1/Rmu = 1/Rt+1/Ry = 6 mS, so Aypi= 
2m/[1/Ritgati{2nf)Ci] and Avpj= 8m/[1/Ritgati(2nf)Cir] are approximately simplified to —Ayp = 


2m/(1/Ri +24) = gmRa, Zur to Rui, and Caj to Ca. Thus, 
a _ i}! (ole m l Ry Een ( l ) Ey ( l as 
T 1)” (0) © (4(2af)Ca) ARa (IOS )Ca lek V2 ex) AR (84.8) 
( 1 eal 1 = 1+ Ra /R; ji 
j(2nf)C,] AR, (i(2af)C, Ar (2m )(Ca + CR /Ra)) Em 


By Wiener-Khintchine theorem, ignoring the small quantities, such as the thermal noise of R1, 
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Ro, and Rr, 


ey =e; +{(R IR) Ea +(e te te; ei) |/4+ (1+ Rin /Ra) ile > (s4.9) 
= Rè e etete oA 
«2 2 2 a 
Ly = P+(2 mf) on L Ta Sn wf) On cy ee mf) [c+ RÓ C Ja 2 5 (s4.10) 


fa Roe, etete R R ) 2 
ei = hina) = = -eje E r E 7 Z [+| C+ E a DE 2 . (s4.11) 
vP d m 


Here, e is the equivalent input noise voltage PSD of the CNRS-HEMT CH and z2 is their 
equivalent input noise current PSD. e? is the equivalent input noise voltage PSD of the Rear-OPA 


and P is its equivalent input noise current PSD. e and e are the thermal noise voltage PSD 


of Rei and Rez respectively. 
e? =e?=4k,T, R, =1.66 (nV)'/Hz , where Tr is 300 K. e?=2.25(nV)?/Hz and 
2 =4(pA)’/Hz in f> 10 kHz [S4R4]. In Eq.(s4.9)-(s4.11), (Rẹ, /R3)eku is two orders of 


magnitude smaller than etete. 1+R,/R,) i2/g2 is one order of magnitude smaller than 
8g a 1 2 HL d a Sin g 


e in Eq.(s4.9). (Ca +C, Ru /Ra y F gi is one order of magnitude smaller than 
Ci (etete; )/& in Eq.(s4.10). (Cy+C, Ra /Ry)(1+ Ru [Ry e/g? is one order of 
magnitude smaller than C, (terte) / A’, in Eq.(s4.11). Further ignoring the small quantities in 


Eq.(s4.9)-(s4.11), 


e+e, +(e te) +65 /4 (s4.12) 
pe as. (2af? C Alen +e +e; e) A ; (s4.13) 
au =-j(2rf)C, (etete) A - (84.14) 


84.2.2 Equivalent input noise voltage and the equivalent input noise current of of TIA 
For TIA, putting Eq. (s4.12) into Eq. (s4.1), putting Eq.(s4.13) into Eq. (s4.2), and putting Eq. 


(s4.14) into Eq. (s4.3), its equivalent input noise voltage PSD r and its equivalent input noise 


current PSD A , its equivalent input noise voltage-current PSD e,i% , and its equivalent input noise 
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current-voltage PSD ipe; are obtained as 


ee, +(eta ter) /A., (84.15) 
ip = 4kpT /Rp + en R? +(er ter ter) (ARRE) tin +(2af) CÈ (etete) /4h , (416) 
ei =(iner) = a/R (Ha A(R) 2) (r) (54.17) 


84.2.3 The equivalent input noise current of the proposed CryoSTM-TIA 
Putting Eq.(s4.15), (s4.16), and (s4.17) into Eq.(s4.4), the equivalent input noise current PSD of 


the proposed CryoSTM-TIA i2 is obtained as 


È =i} +4k,T/R, +(I/R, +1/R,) fei + (cite rei) Ae ar) E +C? (citeres) / 4 | ; 
(s4.18) 


where Cy = Crt+C; and C= Cat+C}+Cj. Eq.(s4.18) is Eq.(3.11) in Article. And, 


2 =i +4k T/R, +(1/R, +R.) eè +i, , (s4.19) 


where pe = (2nf) E e +C G +e +e; è] . Eq.(s4.19) is Eq.(3.12) in Article. 


[S4R1] A. van der Ziel, Noise in Solid State Devices and Circuits, Wiley-Inter-Science, New York, 
(1986). 

[S4R2] Y.X. Liang, Low-noise large-bandwidth transimpedance amplifier for measuring scanning 
tunneling shot noise spectra in cryogenic STM and its applications, Ultramicroscopy, 234 (2022) 
13466. 

[S4R3] Z.H. Qian, Study on noise models algorithms and matrix descriptions for integrated circuits, J. 
Northeast Norm. Univ. 35 (2003) 41. 


[S4R4] Data sheet of THS4021 OPA, https://www.ti.com/lit/ds/symlink/ths4021 .pdf.3. 
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Supplemental file 5: Estimating the DC tunneling current error for CryoSTM-TIA 
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Figure s5-1 For the proposed CryoSTM-TIA, Er = |L-I/|I| vs. R (Re[1 MQ, 1 GQ]) is simulated by TINA-TI with Vi 


=] mV (Black curve) and Vi=5 mV (Red curve). 


For the proposed CryoSTM-TIA, Fig. s5-1 shows the TINA-TI simulation results Er =|Js-|/|J| vs. R 
(Re[1 MQ, 1 GQ]) with Vi=1 mV (Black curve) and Vi=5 mV (Red curve). The simulation results 


show Er<30 ppm, consistent with the calculated results in Article. 
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Supplemental file 6: Estimation of Inv-Amp input offset voltage 


Figure s6-1 DC circuit of the inverting amplifier. 


Table s6-1 Parameters of THS4021 as the Rear-OPA in Post-Amp 


Open loop DC voltage gain aao 97.5 dB 
Input offset voltage drift |Vord/°C 15 pV/°C 
Input bias current Jep, Jon 3 pA 
Input bias current drift |ipp|/°C, |épal/°C Estimated as 100 nA/°C 
CMRR 95 dB 
Ac = Aao/CMRR 2.5 dB 


The Pre-Amp in the left dashed box of Fig.s6-1 is made of CNRS-HEMT CH. The source of CH is 
grounded via a variable resistor Rs. The source of CH is connected to Ri. Ri=R2=1 KQ, Ru =R1, R= 
R2tRp, and variable resistor Rp can be adjustable from 0 to 10 Q. Rui ~ Ri2~Ri_=1 KQ. Its input (the 
gates of CH) is grounded, and its two outputs O; and O2 are suspended. 

Adjust Rp, Rs and Ru, to let CH at the ideal operating point (las, Vas)=(1 mA, 100 mV), and the DC 
voltages at O is 0. The gate-source voltage Ves = —JasRs, and Ve=Vas-Vzs. I assume Vgs =-50 mV for the 
ideal operating point, so Rs=50 Q and Ru= 150 Q. 

For the inverting amplifier, assuming a output voltage drift of vo for 1 °C, the voltage drift at the 
output O; is vı, and it at O; is v2. The voltage drift at the output of the constant-current source is vs, the 
voltage drift at PVR is vp, and the voltage drift at the source of CH is v;. By the following equations, 
the input offset voltage drift of the inverting amplifier can be obtained. 


(vovi) RtH vv /RL= (vpo vi)/ Ri 5 (s6.1) 
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(vt vi1)/RL= ga(v1 Vs) &mVs Ibp > (s6.2) 


(Vo-v2)/Ret (vi-v2)/RL= v2/Ru +ibn , (s6.3) 
vs/Rs = ga(Vi-Vs) —ZmVs 5 (s6.4) 
Aa0(V1tVott-V2) + Ac(VitVott +V2)/2 = Vo , (s6.5) 
Vos = Vo/@ao . (s6.6) 
avo =~ 2mR/(gaRu) . (s6.7) 


where, vos is the input offset voltage drift of InvAmp for 1 °C. From the above equations, vos and vs 
can be solved out as vos ~ ipp(1+gmRs)/gm7lngaRu/2m-Vp(1+emRs)/(2gmRt)-Vor( 1+2gaRit+gmRs)/(22mRr) 
and vs ~ vpRs/(2Rt)~ ibpRs. 

\vp|S<0.4 mV is the current drift of the constant-current source for 1 °C. Hence, |vs|<10 uV, and 
|vos|<17 uV. i.e., in Article, the drift of the gate-source voltage of CH in the inverting amplifier is less 
than 10 uV/°C, and the input offset voltage drift of the inverting amplifier is less than 17 uV/°C. The 
power of PVR is less than 40 mW, and the power of THS4021 is less than 200 mW. It is easy to control 
the temperature fluctuations of PVR and that of THS4021 within 0.01 °C with the temperature Control 


System, so the total input voltage drift can be controlled within 170 nV. 
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